Determination of the ttE scattering lengths from the weak decays of Ac 
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The scattering lengths of the ttS systems are key quantities in understanding the structure of 
the A(1405) resonance and the subthreshold extrapolation of the KN interaction. We demonstrate 
that the ttE scattering lengths can be extracted from the threshold cusp phenomena in the weak 
Ac — ttttE decays, analogously with Cabibbo's method for determination of the tth scattering 
length. We show that the substantial cusp effect should be observed in the spectrum, when the ttE 
interaction in 7 = is strongly attractive to generate a near-threshold singularity, such as a bound 
state or a virtual state. 
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I. INTRODUCTION 

The possibility of kaon bound states in nuclei [H 
has been intensely discussed recently [3l-[l0j and is one of 
the central subjects in forthcoming experiments at Japan 
Proton Accelerator Research Complex E15 FOPI at 
GSI [12, and AMADEUS at DA$NE An impor- 

tant theoretical concept to describe kaonic nuclei is the 
two-body KN-nT, interaction, which closely relates to 
the description of the A(1405) resonance in the 1 = 
amplitude. Theoretical models of the KN-ttI] scatter- 
ing amplitudes have been constrained by the total cross 
sections of low-energy K~p scattering into various final 
states and the threshold branching ratios of the K~p 
channel. Recently, due to intensive experimental activ- 
ities, new data of the ttE invariant mass distributions 
are becoming available, and the analysis of the energy 
level and width of the kaonic hydrogen provides pre cise 
information regarding the K~p scattering length [IJl . 

In this way, there are many experimental data around 
the KN threshold, while the amplitude at far below 
the threshold is not well constrained. This means that 
the KN-nYi interaction relevant to the study of the K- 
nucleus systems with strong binding is achieved only 
through the subthreshold extrapolation from the KN 
threshold. For instance, the chiral unitary model (l5l - [l9l | , 
a coupled-channels approach with the low-energy inter- 
action constrained by chiral symmetry, has revealed an 
interesting two-pole structure [l^. This indicates that 
the nominal A(1405) resonance is not a single state but 
a superposition of two states with different properties. It 
is shown that this two-pole structure has substantial in- 
fluence for the KN interaction [U and the KNN-ttY.N 
system (T^]. However, it is found that the position of 
the lower energy pole is sensitive to the details of the 
model, due to the lack of the information of the ttS in- 
teraction [H^. 

The scattering length characterizes the property of the 
two-body interaction. The attractive/repulsive nature of 
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the interaction is reflected in the sign of the scattering 
length. The magnitude is sensitive to the singularity of 
the amplitude close to the threshold. If there is a bound 
(virtual) state, the scattering length becomes large and 
negative (positive), as is known for the ^Si Q'Sq) NN 
scattering [2^. In the case of the KN-ttT, interaction, 
it is shown that the position of the lower energy pole 
of A(1405) is closely related to the ttE scattering length 
with isospin 7 = [23|. Therefore, it is highly desirable to 
determine the 7rS(7 = 0) scattering length to impose fur- 
ther constraints on the low-energy KN-irYi interaction. 
There arc several ways to extract the hadron scattering 
length in experiments. The standard method is to mea- 
sure the energy shift and the width of the hadronic atom 
state [13, [Ik l2IJ 01" to extrapolate the low-energy phase- 
shift data down to the threshold [12 ■ Unfortunately, 
both methods cannot be applied to the ttS case. 

Recently, an alternative method to determine hadron 
scattering lengths was proposed by Cabibbo [H, H^] , in 
which the cusp effect at the tt+tt" threshold is used to 
extract the tttt scattering length in the K^ — >■ Tr+Tr^Tr" 
decay. The structure around the cusp at the tt+tt^ 
threshold is related to the scattering amplitude of the 
TT^^TT^ — > tt'^tt'^ process, which is nothing but the tttt 
scattering length (the relationship between the scattering 
length and the cusp phenomena in the tttt spectrum has 
been discussed also in Ref . [s^ ) . The cusp effect is indeed 
observed in the high-statistics experimental data [1^ Is^l 
which leads to the determination of the tttt scattering 
length. 

In this paper, we propose to apply a similar methodol- 
ogy to the ttS scattering lengths by use of the Ac — )• ttttS 
decays. This decay process has been measured as the 
Ac peak in the three-body invariant mass spectrum in 
Refs. [H-[42. The branching ratios to the relevant decay 
channels are reported as a few percents by the Particle 
Data Group (PDG) [12 , so the extraction would be fea- 
sible at Belle and Babar where a huge amount of Ac is 
produced. 

The contents of this paper are as follows. In Sec. |lll 
we discuss the general property of the threshold cusp 
phenomena in the Ac ttttS process. We introduce 
an expansion scheme to extract the scattering length. 
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keeping the value of the scattering length unspecified. In 
Sec. Illli we present theoretical estimations of the decay 
process and show the expected spectra for several values 
of the scattering lengths. The last section is devoted to 
the summary. 

II. ANALYSIS OF THE CUSP PHENOMENA 

Here we analyze the threshold cusp phenomena in the 
Ac — > ttttI] decay. We review the general argument of the 
determination of the scattering length from the thresh- 
old cusp effect and summarize the possible ttE channels 
in the A^ decay where we have a chance to observe the 
cusp structure. We then discuss the method to extract 
the scattering length with minimal model assumptions 
by expanding the amplitude in terms of the momentum 
variable. We also present an extension of the framework 
to the case with complex amplitudes. 



which is depicted in Fig. [Ijb) . Because of the mass differ- 
ence, the (ttE)^ threshold appears in the mass spectrum 
of the (ttS); channel slightly above the (ttS); thresh- 
old. The cusp phenomenon occurs at this higher energy 
threshold. 

It is important to note that the amplitude in Fig. [Ub) 
contains the vertex of (ttS)/! (ttE);. The threshold 
cusp takes place with the vanishing momentum of the 
(ttT,)^ state, and it contains the information of the on- 
shell scattering amplitude of the (7rE)/i — >• (ttE)/ process. 
Following Ref. [U, we define the off-diagonal ttS scat- 
tering length as (the real part of) the amplitude fh^i at 
the higher energy threshold W = Wth'- 

ah^i=fh^iiW ^Wtb), (1) 

which is reflected in the behavior of the cusp structure. 
In the following, we show how to extract the scattering 
length ll]) along the same line with Ref. jssj . 



A. Ac ttttE decay and the ttE scattering length 



B. Possible decay modes for the scattering length 



Let us consider the weak decay of Ac into two pions and 
one S baryon. We denote a ttE pair in the flnal state as 
(ttE); and concentrate on the threshold energy region of 
this pair. In this case, the additional primary pion has 
a large momentum because of the large phase space of 
the decay process, as we see in Appendix [X] Choosing 
an appropriate charged state for the (ttE); pair, other 
charge combinations of the ttE state may be allowed at 
slightly higher energies than the (ttE); threshold due to 
the isospin violation in the particle masses. We denote 
the higher energy channel of the ttE pair as (7rl])/j. The 
specific charge states of (ttE); and (7rE)/j will be given 
in Sec. Ill B1 in connection with possible experimental ob- 
servations. 

The dominant part of the Ac decay is given by the 
direct process: 

Ac ^7r(7rE)i, 

as shown in Fig.[lja). In addition, we may have the final- 
state interaction term with the intermediate (ttE)/! state 
as 

Ac ->7r(7rl])ft 7r(7rE)i, 



(a) 



(b) 



FIG. 1. Decay diagrams for the Ac — i- 7r(7rS); process: (a) 
direct decay and (b) decay through the (ttE);, state. Solid 
circles denote the weak process of Ac — ttttE and the open 
circle represents the scattering amplitude of the (ttE);,. — >■ 
(ttE); process. 



We first consider the energy levels of the ttS channels 
with different charge combinations to specify the (ttS); 
and (ttE)/! channels. Since the charge Q = ±2 states 
have no coupled channels, they are not relevant in the 
present discussion. The threshold energy levels for the 
Q = 0, -1-1, and —1 channels are shown in Fig. [21 Unlike 
the TTTT case where the masses of tt^ and tt^ are identical, 
E~ is heavier than E"*", so the ttE channels show a rich 
spectrum. Among them, there are about 10-MeV mass 
differences in the following transitions: 



tt^EO, 



^+E0 ^^"E+ 



1335 



1330 



1325 



1320 



A- 



7t"Z" 



Q=0 



Q=+l 



Q=-l 



FIG. 2. (Color online) Threshold energies of the ttE states 
with charge Q. The arrows indicate the possible decay modes 
in which visible cusp structure may appear. 



3 



TABLE I. Decay branching ratios Ti/T from PDG [43 



mode 


Ac 7r(7rE)h 


Ac 7r(7rE), 


a-+ 


1.7 ±0.5 % 


3.6 ± 1.0 % 


a°o 


1.7 ±0.5 % 


1.8 ±0.8 % 


a»+ 


1.8 ±0.8 % 


Not known 



percents. If the precise measurement of the ttE spectrum 
around the threshold is performed, it will be feasible to 
extract the scattering lengths by the method explained 
in Sec. HTD] 



C. Threshold cusp effect 



We expect to observe the cusp effects by regarding these 
decays as {TrT,)h — > (ttS);. 

Using the isospin decomposition with the phase con- 
vention given in Eq. (|B6p . we can express the scattering 
lengths in these channels as 

=ia" - ifl^ + H , 



,00 



where the scattering lengths arc labeled by the charges of 
the final states, is the scattering length with isospin / 
in the isospin basis, and the ellipses represent the isospin 
breaking corrections^ Although we have three equations, 
they are not linearly independent if we neglect the isospin 
breaking corrections as 



,00 



-,0+ 



-,00 



and aP^ 



Therefore, even if we extract a" 

periments, in order to determine all three a^s, we need 
one additional input. We may adopt chiral perturbation 
theory (ChPT) or the lattice QCD simulation for the de- 
termination of a^^^ where the interaction is presumably 
weak and repulsive. 

Let us check the decay modes of Ac which lead to these 
rescattering amplitudes. Since the charge of Ac is Q = 
+1, the possible decay modes are 

:Ac^7r+(7r+S")^7r+(7r-I]+), (2) 

a°":Ac^7r+(7r+E-)^7r+(7r"S"), (3) 
o+.A ^ ^o/_+vOn ^ Jv+^ 



a 



(4) 



For later convenience, we refer to these modes as the 
extracted scattering lengths, i.e., Ac — > 7r+(7r+I]^) — > 
7r+(7r~E+) is called mode a '". The experimental 
data [i^ of the branching ratios r,;/r for these channels 
are summarized in Table HI 

The mass spectrum is extracted from the Ac — > 7r(7rl])( 
mode, while the Ac — > 7r(7rS)/i mode is used to normal- 
ize the strength of the amplitude, as we will see below. 
Fortunately, most decay modes are experimentally ob- 
served, with the branching ratio of the order of several 



^ In principle, can be complex, since the ttA channel is open 
at the ttS threshold. However, the transition 7rS(/ = 1) — > ttA 
vanishes at the leading order in chiral perturbation theory. In 
the following we assume the imaginary part is negligible. 



Let us consider the Ac decay process in Fig. [1] paying 
attention to the imaginary part of the loop function. To 
concentrate on the cusp phenomena, in this section, we 
simply ignore the diagrams other than those in Fig. [TJ 
assuming that the other final-state interactions are slowly 
varying in the relevant energy region. The decay width 
of the process Ac — )■ 7r(7rS])/ can be calculated as 



where M is the relativistic scattering amplitude of the 
process, dUs is the three-body phase space and E denotes 
the spin summation. We denote the mass of the E^ (tt/i) 
by Mh (to/i), and the invariant mass of the (ttE)/ system 
by W. The mass spectrum of the decay with respect to 
W is given by 



dT 
dW 



dfl E|7W|^ 



with dn ~ dli^/dW. Since we are interested in the cusp 
structure which arises from the nonanalytic behavior of 
W, we concentrate on the amplitude as a function of W, 
assuming that the dependence on the other kinematical 
variables can be factorized. Integrating the other vari- 
ables, we arrive at the expression 



dT 



M,W 



dW Idn^MA^ ' 



E|X|2 / duJie{l^A^), 



(5) 



A 



Pi 



Ml 



^2 



Pi = 



where Ma^ is the mass of the Ac, (wtt) and pi (p^) are 
the energy and momentum of tt; (the primary pion), and 
Q{x) is the step function. 

To appreciate the threshold cusp effect, we focus on 
the property of the loop function G{W) in Fig.[TJb). The 
loop function can be written in the spectral representa- 
tion as 



G{W) 



1 

'2^ 



w 



where the threshold energy is Wth = 
phase-space factor p{W) is given by 



± (subtractions), 
Mfi + nih and the 



piW) =2Mh 



qjW) 
inW' 
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with the three-momentum function 



q{W) 



2W 



The real part of the loop function G{W) depends on the 
subtractions, but the imaginary part can be determined 
only by the kinematics. Since the scattering amplitude 
with the diagram in Fig. [T] contains the loop function, we 
obtain the condition for the amplitude 



Im7W(Ty)aIm G{W) 



Q[W-Wtx,). 



(6) 



This means that the imaginary part of the amplitude 
appears suddenly at the threshold. 

Considering the two diagrams shown in Fig. [l] we can 
decompose the Ac — >■ 7r(7rE)i amplitude into two parts: 



M{W) =Mo{W) + iMi{W)q{W) for W > Wth, 



(7) 



with Mo{W) and Mi{W) being analytic functions of W 
at the threshold Wth. Here we also assume that Aio{W) 
and Aii{W) have no imaginary part, and the extension 
to the complex amplitudes is discussed in Sec. Ill El The 
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amplitude ^ satisfies the condition ([6]), because the 
(ttS)/! loop is the only source of the imaginary part in 
this process. Physically, the amplitude of the direct pro- 
cess (a) and the real part of the indirect process (b) are 
included in A^O; a-nd the amplitude corresponding to the 
imaginary part of the indirect process (b) is included in 
iMiq{W). The function q{W) can be analytically con- 
tinued to < Wth where it becomes pure imaginary. 
For later convenience, we define the dimensionless quan- 
tity 



th 



W^) iq{W) 



(8) 



so 5 is real (imaginary) below (above) the threshold and 
,52 = _|^|2 foj. 1^ > W^thS The amplitude then can be 
rewritten as 



M{W) =Ma[W) +Mi{W)mh5. 



(9) 



Because of the property of (5, the second term of Eq. © 
interferes with the first term below the threshold, while 
such an interference does not occur above the threshold. 

Using the expression © and M* = Mq + MimnS*, 
we calculate the amplitude square as 



(Mo)^ + iMimn)^ for W > Wth 

{Mof + 2MnMimhS + {Mimh)^S^ for W < Wth 



(10) 



r 



Note that 6 vanishes at the threshold, so the spectrum is 
continuous a.t W = Wth- On the other hand, its deriva- 
tive with respect to W is not continuous; 



due to the kinematical condition, but the simple inter- 
ference pattern in Eq. (|T0|) is realized only for the real 
valued Mo amplitude. 



d\M\^^ 



dW 



d\M\^^ 



oc 



2MoMimhMh 1 
Mh + rrih 5 



dW 



0{5). 



This means that the (ttE); spectrum is continuous but 
not smooth at the (ttS)/! threshold. This singular be- 
havior is called the threshold cusp. Although we have 
assumed that AAq is real in the present case, the mech- 
anism of the threshold cusp is the same even if A^o has 
an imaginary part, since the discontinuity is caused by 
the nonanalytic term 5 as a function of W , as shown in 
Sec. Ill El Thus, the cusp always occurs at the threshold 



^ There is another branch point at ly = My^ — , but we will not 
consider this because the energy region cannot be reached in our 
analysis of vrS spectrum around the threshold. 



D. Expansion of the amplitude and spectrum 

In order to extract the scattering length, it is desirable 
to parametrize the amplitudes around the threshold Wth- 
The aim of this section is to relate the ttE scattering 
length to the near-threshold behavior of the spectrum. 
Since we concentrate on the energy region around the 
(ttS)/! threshold, the variable 5 should be small. As we 
will see below, the relevant energy region is 1 14^ — Wth\ ^ 
10 MeV, where the value of 5 is \5\ < 0.37. This means 
that the 0{S^) error is about 5 % level. 

Since the Mo amplitude is an analytic function of W, 
it has only even powers in the S expansion: 



Mo{S) =A^r 



(11) 



The Ml part of the amplitude is related to the scatter- 
ing length. Identifying this contribution as the imaginary 
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part of the diagram in Fig. [Ijb) , we can write the amph- 
tude above the threshold as 

where Aig is the amphtude of the weak Ac — > tt{ttY,)ji 
decay. The amphtude Aih^i represents the transition 
amphtude for the (ttE)/! —> (ttS); process and fh^i is 
its nonrelativistic counterpart. Note that aU the final- 
state interactions via the strong interaction are implic- 
itly included in fh^i- The Mq amplitude is an analytic 
function of W so 



(12) 



The low-energy behavior of the nonrelativistic scattering 
amplitude can be parametrized by the scattering length 
ah-^l■ According to Eq. ([T|), the scattering length is de- 
fined at the higher energy threshold, so we expand it in 
terms of 5 as 



1 



ah^i 



- iq + ' 



^ah->i + al^iTUhS H . (13) 

Using Eqs. (fT2l) and (fT3| . we can express Mi a.s 

Mi{S) ^Mo'-"^ah^i + Mo''''^m,,al^i6 + 0{5^). (14) 

Combining Eqs. ©, pT|) . and (HH), the expansion of M 
up to is given by 

M{S) =mIj"^ +Mo''"'>mhah-,iS 



The amplitude square is then given by 



A + C'\S\^ + Oi\S\^) 
A + BS + CS^ + 0{S^) 



for W > Wth 
for W <Wth' 

(15) 



where the coefficients are given by 



h^l 



These coefficients can be extracted using the experimen- 
tal spectrum, the three-body phase-space factor, and the 
fitting by a polynomial of i5 around the threshold Wth- 
Since the B coefficient is proportional to the scattering 
length, we can extract its absolute value as 



\ah^i\ 



\B\ 



2mnVA\M"a 



,'i(0)| 



and the sign is given by 



\ah^i\ 



'\mT\ 



B 

m 



The coefficient can be obtained as the leading 

contribution for the Ac — > 7r(7rS])/j process around the 
threshold. It is also possible to calculate TWq when the 
information of the weak process is well under control. If 



we determine the relative sign of A^q"'' and A^q*^'^'', the 
sign of the scattering length a^^i is determined by the 
sign of the coefficient B. This information, together with 
the fitting of the (ttE); spectrum in the Ac — > 7r(7rl]); pro- 
cess below and above the (ttE)/! threshold by Eq. ([T5|), 
leads to the determination of the scattering length at^i- 



MO) 



E. Extension to the complex amplitude 



So far we have assumed that the weak interaction am- 
plitudes and Ml are real. This is valid at the lead- 
ing order, but several final-state interactions may bring 
an imaginary part as we discuss in Appendix [Xj Here we 
consider the possible modification of the formulation by 
the imaginary part of these amplitudes. 

Since the overall phase of the amplitude does not con- 
tribute to the spectrum, we consider a relative phase 
between Mq and Mi amplitudes. The imaginary part 
other than Eq. (jH]) is a smooth function of W , so we in- 
troduce a constant phase e** to generalize Eq. © as 



M{W) =Ma{W) + Mi{W)e'%%h8. 



with Mo and Mi being real. In this case, the amplitude 
square is given by 
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(Mo)^ + 2MoMimh\S\sin0 + {Mimh)^\S\'^ for W > Wti, 
(Mo)^ + 2MoMimhScose + {Mimh)^S^ for W<Wtb' 



In this case, the Uncar \S\ term remains in the spectrum 
above threshold. 

The expansion of the amphtude square in terms of |5| 
is given by 



\M\' = 



A + B'\S\+C'\6\^ + 0{\S\*) for W > Wt 



th 



A + BS + CS^ + 0{6^) 
where the coefficients are given by 

B ^2M^o^Mo''°^mhah^i cos0, 



for W <Wth' 



B' ^2Ml^^M'^^°^mhah^i sin 



61 = B tan 6*, 



C =2Mi''^Ml^^ + 2M[l'^M'i^"^mlal_, cosO 
C'=- 2M^°^M'i^ - 2M'i^MT'^mlal^i coat 

The scattering length is 

|B/ cos 6*1 



\ah~^i\ 



2mhVA\Mo 



M 



(0) 



M 



B/cost 



\mT''\ W/cosey 

Thus, the scattering length can be extracted in the same 
way as before when the relative phase between M.q and 
Ail is under control. In practice, the relative phase 
mainly stems from the final-state interactions of the tttt 
scattering, which are well understood by the phase-shift 
analysis. 

Even if the relative phase is not known in advance, the 
magnitude of the scattering length can be extracted from 
the coefficients of the linear terms: 



^B^ + {B'Y 
'2mh^/A\MT\' 



\ah^i\ 



In addition, cos 9 can be determined from the 5^ terms if 

(2) 



we neglect the term 



cosfc'= — — for 



(2) _ 



In this way, extraction of the scattering length is possi- 
ble for the complex amplitudes, when information of the 
final-state interaction and/or high-precision data of the 
decay spectrum are available. 



III. ESTIMATION OF THE MASS 
DISTRIBUTION 

So far we have discussed how to extract the ttS scat- 
tering lengths from the observed decay spectrum. In this 
section, we present theoretical estimation of the mass 
spectrum (|15p , in order to examine the sensitivity of the 
cusp effect on the ttE scattering lengths. 



A. Input parameters 

To evaluate the mass spectrum, we need to determine 
the amplitudes for the weak process and the value of 
the scattering length. We first consider the weak pro- 
cess of the Ac — >■ TTTrE. Since we are not aiming at the 
construction of the realistic decay spectrum, we simply 
approximate the amplitude by the leading-order constant 

Mq = A^o"^ and = A^o'°\ which would be sufficient 
for the present purpose. We first put 9 = Q and the ef- 
fect of the relative phase will be considered in the end 
of this section. The magnitude of these amplitudes are 
determined by the partial decay widths derived from the 
central values of the lifetime of the Ac and the branch- 
ing ratios given in Table IH The results are summarized 
in Table |lll For the mode, the branching ratio to 
the lower energy channel is not known and we assume 
the same strength with the higher energy channel. In 
addition, the relative signs of the amplitudes have been 
determined by evaluating the quark diagrams of the weak 
decay in Appendix |B] In the following, we adopt the cen- 
tral values for the numerical calculation of the spectra. 

Next we consider the ttS scattering lengths. It is 
known that the scattering length is well described 
by the leading-order term in ChPT, up to 10%. Thus, 
as a first trial, we may adopt the leading-order ChPT 
to estimate the ttS scattering lengths. Higher-order con- 
tributions to the scattering lengths are calculated in the 
heavy baryon formalism [4J] and in the covariant formal- 
ism with the infrared regularization scheme (45| . With 
the definition ((TJ, the leading-order term in ChPT gives 



TABLE II. Strengtlis of the weak interaction vertices. 



mode 


TWS (10^^ MeV~') 


Mo (10-^ MeV-i) 


a- + 


1.46 ±0.43 


2.10 ±0.58 




1.46 ±0.43 


-1.49 ±0.66 




-1.49 ±0.66 


-1.49 ±0.66^ 



^ Assumed to be the same with JVIq 
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the ttE scattering length for the channel /i — ?> / as 



4/2 



2{Mh+mh) 

where / = 92.4 MeV is the pion decay constant and 
the ellipsis denotes the isospin breaking correction. The 
coefBcicnt Cu^i is the group theoretical factor [4^, |43] , 
which takes the values 0, 2, and —2 for the a '", a™, and 
aP^ modes, respectively. Thus, the leading-order ChPT 
predicts 



=0.00 fm, 



,00 



0.23 fm. 



-0.23 fm. 



In the isospin basis, the coupling strengths are given by 
4, 2, and —2 for / = 0, 1, and 2 channels. The strength 
C = 2 is the same with the TrN{I = 1/2) channel, so 
we expect that the ChPT prediction works well. This is 
indeed the case for the 1 = 2 scattering length, since the 
lattice QCD simulation by NPLQCD Collaboration ^ 
found that a7r+i;+ = —0.197 ±0.017 fm, which is compa- 
rable with the leading-order ChPT prediction of —0.2294 
fm. On the other hand, C = 4 in the / = channel is 
twice as strong as the ttN interaction so the rcsumma- 
tion effect may be important. Indeed, it is found that the 
chiral 7rS(/ = 0) interaction generates a resonance state 
by the single-channel resummation [2l[ . The relation be- 
tween the pole singularity around the ttS threshold and 
the ttS scattering length is studied in Ref. [l^l, which 
shows that the position and the nature of the pole are 
sensitive to the value of the scattering length. Depend- 
ing on the nature of the lower pole of A(1405), the typical 
ttY,{I = 0) scattering lengths are found as 



a^-° = i 5.5 fm 



1.2 fm: resonance with E = +31 McV 

virtual state with E = —6 MeV 
-2.3 fm: bound state with E = -10 MeV 



where E denotes (the real part of) the pole singularity 
measured from the threshold[f| It can be seen that the 



TABLE III. Estimation of the scattering lengths in units of 
fm. ChPT stands for the leading-order prediction in chiral 
perturbation theory. The other results are obtained in the 
resummation scheme for / = amplitude when the pole sin- 
gularities are developed at the energies in parenthesis away 
from the threshold [2711. 







a- + 




a"+ 


ChPT 


0.46 


0.00 


0.23 - 


0.23 


Resonance (-|-31 MeV) 


1.20 


0.25 


0.48 - 


0.23 


Virtual state (-6 MeV) 


5.50 


1.68 


1.91 - 


0.23 


Bound state (-10 MeV) 


-2.30 - 


-0.92 


-0.69 - 


0.23 


3 The calculation in Ref. [27 


was done 


in the 


isospin limit. 


The 



resonance case is estimated by taking an average of two different 
solutions, Al and B E-dep in Ref. [27| . 



(a) a * mode 

a - +0.25 fm (Resonance) 

a = 0.00 fm (ChPT) 



(b) a ^ mode 
- - a = +1.68 fm (Virwal) 

a = -0.92 fm (Bound) 

a = 0.00 fm 




1330 1340 
W [MeV] 



1330 1340 
W [MeV] 



FIG. 3. (Color online) Spectra of the 7r~E^ channel in the 
Ac — > 7r"'"(7r~E"'") decay with several values of the a""*" scat- 
tering length. 



scattering length takes a large value when the pole is 
located close to the threshold, and the sign of the scat- 
tering length is inverted when a bound state appears. 
With these values for the / = scattering length, while 
the 1 = 1 and the 1 = 2 components are fixed as the 
values by the leading-order ChPT, we obtain the scat- 
tering lengths by use of Eqs. ([5]), ([3]), and @ as shown 



in Table IIIIl Even when 



,/=o 



is as large as 5.50 fm, the 



scattering lengths a and a'^° are not much enhanced, 
because of the 1/3 factor in Eqs. ([2|) and ([3]). Note that 
a°+ does not depend on the / = amplitude as observed 
in its isospin decomposition 



B. Examples of ttS spectrum 

Using the information given in the previous section, we 
calculate the (ttS); mass spectrum in the Ac — 7r(7rE)/ 
process by use of Eq. ([5]) . In Fig. [3l we present the mass 



spectra of tt E+ channel in the A^, 



r+(7r E+) mode 



for several values of a . The dotted lines represent the 
case when the scattering length is set to be zero. The 
spectrum shows the pure phase-space distribution and 
no cusp structure appears, because of the absence of the 
B term in Eq. In this mode, a = 0.00 fm also 

corresponds to the prediction of the leading-order ChPT, 
as shown in Table IIIIl 

When a finite scattering length is adopted, the cusp 
structure appears (the solid lines and the dashed line). 
The cusp effect is not very significant for the case with 
a ^ = 0.25 fm. but clear structures are visible for a ^ = 
1.68 fm and a ^ = —0.92 fm. In addition, the sign of the 
scattering length is reflected in the behavior of the lower 
part of the spectrum in comparison with the phase-space 
distribution. 

Next, in Fig. [H we show the spectra of the Ac ^■ 
7r+(7r°S°) decay and the Ac n°{TT°^+) decay with a°° 
and a^"*", respectively. Again, the dotted lines represent 
the case when the scattering length is set to be zero. We 
observe that the prominent cusp structure is seen when 
the scattering length is as large as 0.5 fm. Because of 
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FIG. 4. (Color online) Spectra of the tt^E" channel in the Ac — > tv'^{tv'^TP) decay [(a) and (b)] and the tt^E"*" channel in the 
Ac — )■ 7r''(7r''E+) decay (c), with several values of the scattering length. 



the difference of the relative sign of the Aio and Mq 
amplitudes, the interference pattern in the a*^*^ mode is 
inverted from that in the a ^ mode (see Table HI)) . 

Let us compare three decay modes. In both the a 
and a^" modes, the effect of the / = component is re- 
flected in the cusp structure, so these modes can be used 
to extract the isosinglet scattering length. A nice feature 
of the a ^ mode is the large branching fraction of the 
decay (^3.6%) which is advantageous for the experimen- 
tal observation. On the other hand, the scattering length 
vanishes in the leading-order ChPT, and, hence, the ex- 
pected scattering length a ^ by the resummation effect 
is not very large. The a^" mode has a smaller branching 
ratio to the final state, but the value of is in general 
expected to be larger than a as shown in Table IIII[ 
which causes a stronger cusp effect. The a^~^ mode is 
not useful to extract the 1 = component, and the ex- 
pected cusp structure seems to be relatively weak. In 
addition, observing two neutral pions in the final state 
would be experimentally difficult. Nevertheless, once the 
precise decay spectrum is observed, this mode serves as 
a precise test of the ChPT prediction for the scattering 
lengths in / = 1 and 1 = 2 channels. 

In order to estimate the accuracy of the spectrum mea- 
surement for the feasible determination of the scattering 
length, we consider the deviation of the spectrum from 
the phase-space distribution 



dev(a) 



dT/dW{a) ~ dT/dW{a = 0) 
dT/dW{a = 0) ■ 



(16) 



For a given a, the magnitude of this deviation provides 
the lower limit of the experimental accuracy to extract 
the scattering length. We also note that the deviation is 
always zero at the higher energy threshold Wth = + 
TO/i, so the phase-space distribution can be normalized 
at this point. In Fig. [SJ we show the deviation plot for 
the Ac — >■ 7r+(7r^S+) decay with several values of the 
scattering length a. We observe that accuracy of about 
20% (10%) is required to determine the scattering length 
of \a\ = 1 fm (0.5 fm). 



a mode (deviation) - 

— a = +1.00fm 
a = +0.50 fm 

- - a = -0.50 fm 
.a = -1.00fm 




1330 1340 
W [MeV] 



1350 



FIG. 5. (Color online) Deviations from the phase-space dis- 
tribution defined in Eq. (jlGp for the Ac — 7r"'"(7r~E"'") decay 
with several values of the scattering length a. 



We also study the effect of the relative phase be- 
tween Mq and Ml. As discussed in Appendix 1X1 the 
main source of the contamination process to the Ac — > 



■(tt+E- 



'E+) decay is the tttt rescattering. If 



the primary 7r+ interacts with the n'^ in the intermediate 
state, the Ai^ has an imaginary part. This is the tttt scat- 
tering in the isospin 1 = 2 channel. The relevant energy 
region of the tttt invariant mass is about 400-800 MeV, 
and the phase shift of this tttt channel ranges between 
—5° to —20°, with almost linear energy dependence jioj . 
Although the phase depends on the invariant mass of the 
TTTT pair, we may represent the tttt rescattering effect by 
a single energy-independent relative phase. Thus, for a 
rough estimation, we introduce an energy-independent 
relative phase 9 = —12° and calculate the spectrum with 
the formulas in Sec. Ill El The results are shown in Fig.|6l 
One observes that the spectrum is slightly modified from 
Fig. [3l especially above the threshold. The existence of 
the relative phase affects the interference pattern of the 
two amplitude, while the cusp structure itself remains, 
as discussed in Sec. Ill El 
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(a) a * mode 

a - +0.25 fm (Resonance) 

a = 0.00 fm (ChPT) 



1330 1340 
W [MeV] 



(b) a * mode 
- - a = +l.68 fm (Virtual) 

a = -0.92 fm (Bound) 

a = 0.00 fm 




1330 1340 
W [MeV] 



FIG. 6. (Color online) Spectra of the 7r~E^ channel in the 
Ac — >■ 7r"'"(7r~E''") decay with several values of the a""*" scat- 
tering length with the relative phase 9 = 12°. 
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1.2 


-(a) 


1 1 1 


1.0 
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0.8 
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0.2 






0.0 




1 1 1 



[GeV'] 



IV. SUMMARY 

We study the ttS mass distribution in the Ac — > ttttE 
decays to extract the ttE scattering lengths. The thresh- 
old cusp effect is discussed in detail, and the relation be- 
tween the ttS scattering length and the mass distribution 
is derived. We show that the ttS scattering length can 
be extracted by the expansion of the observed spectrum 
around the cusp structure. 

We then evaluate the ttS spectrum for several values 
of the scattering lengths obtained by chiral perturbation 
theory and the resummation technique. It is found that, 
if the scattering length is as large as 0.5 fm, a substantial 
cusp structure will be observed with about 10 % devia- 
tion in the experimental data. This value of the scatter- 
ing length corresponds to the case when the / = ttE 
scattering amplitude has a bound or virtual state about 
10 MeV below the threshold. Thus, the analysis of the 
cusp structure will give important information on the ttS 
interaction around the threshold. 

Note, however, that we have presented a rough esti- 
mation of the spectrum and there are several other con- 
tributions to the same process which are not included in 
the present calculation. A detailed Dalitz analysis of the 
Ac — ?> ttttS process is highly desired for the calibration of 
the theoretical model of the weak process as well as for 
the investigation of the strong final-state interactions in 
the ttS channel. 

In the Ac decay process, it is possible to provide two 
constraints on the three different isospin component of 
the ttS scattering lengths. To complete the determina- 
tion of all three isospin components, for instance, we need 
to determine the 1 = 2 component by the lattice QCD 
simulation Hll^,!!!. 
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Appendix A: KINEMATICS OF THE Ac ^ ttttE 
PROCESS 

Here we summarize the kinematics of the Ac — tttt/S/ 
process for the assessment of the possible contamination 
processes. We consider the threshold energy region of 
the 7r;S; pair. In this appendix, we adopt the isospin 
averaged masses for S and tt: Ms = 1193 MeV and 
= 138 MeV. 

Let us, first, estimate the typical momenta of the three 
particles. Consider the extreme case when the relative 
momentum of the 7r;S; system is zero and the invariant 
mass of the 7r;S; pair is given by + M^. In this case. 
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the magnitude of the three-momentum of the primary tt 
is calculated as 



(a) 



(b) 



Kl 



747 MeV. 



This is also the magnitude of the total momentum of the 
7r;E/ pair which goes in the opposite direction to the tt 
momentum. With the condition of the vanishing relative 
momentum, the momenta of tt; and E; is given as 



m. 



KI-77 MeV, 



Me + 



\p^ \ - 670 MeV, 



where a large fraction of the momentum of the tt;!]/ sys- 
tem is carried by the E; baryon, because of its heavy 
mass. 

We next show the Dalitz plots of the ttttiJ^i dc- 

cay in Fig. El [Fig. [Ha): M^^, vs. M^^^; Fig.E^b): M^^^ 
vs. M^^^J- The relevant kinematical region for the ttjE; 
threshold production is the left edge of the Dalitz plots in 
these figures. For reference, we put vertical dashed lines 
at the energy (Ms +111^, + 20 MeV)^ Figure [Tjb) shows 
that the invariant mass of the "wrong" ttE; pair for this 
energy region is above 2 GeV, while Fig. ^a) indicates 
that the invariant mass of the tttt; pair stays around 600 
MeV. This difference is caused by the different momenta 
of TT; and E; discussed above. Since E; goes with larger 
momentum than tti in the opposite direction to the pri- 
mary TT, the relative momentum of the E;7r pair is larger 
than that of the tttt; pair. As a consequence, the invari- 
ant mass of the ttE; system is very large, while that for 
the TTTT; system is moderate. 

Based on the kinematical consideration, we discuss the 
possible contamination processes. In addition to the di- 
agrams in Fig. [1] the rescattering of the other combina- 
tion of the final states may occur, as shown in Fig. [5] In 
Fig.m^a) the pions are rescattered, while in Fig.|8l^b), the 
E is rescattered by the primary pion which is not used 
to form the mass spectrum. Since the relevant invariant 
mass of the ttE; pair is above 2 GeV where most of the 
prominent hypcron resonances do not contribute, we may 
ignore the final-state interaction effect of Fig. |SlJb) type. 
In particular, for the a mode, this corresponds to the 
/ = 2 7r+E+ rescattering, which should be small. On the 
other hand, the final-state interaction of Fig. IH^a) may 
be important, since the relevant energy region is near the 
p and a resonances. Because of the isospin, the p (cr) me- 
son can contribute to a ^ and a^" modes (a and a*^^ 
modes) . The a meson is very broad so the energy depen- 
dence of the amplitude is considered to be small. The 
p meson may cause a moderate energy dependence, but 
the upper limit of the branching ratio of the Ac — E+p*' 
is given by < 1.4% |43|. Therefore, one may also as- 
sume that the effect of the p contribution to the present 
argument is minor. 

We also notice that the diagram in Fig. [5] (a) has an 
imaginary part. If we replace the final state TrjE; by 



FIG. 8. Decay diagrams for the Ac 7r7r;E; process: (a) tttt; 
rescattering and (b) ttE; rescattering. Solid circles denote the 
weak process of Ac — ttttE and the open circles represent the 
rescattering amplitudes. 



TTjiE/i, it can be considered as a complex Mi amplitude 
discussed in Sec. Ill El Since the rescattering takes place 
by the primary tt and tt/, the charge combinations are 
given by 



,00 



a°+ : 7r%+. 



In the a mode, the rescattering is in / = 1, which is 
found to be small as < 1.4% [4^. In the a ^ mode and 
the a™ mode, the rescattering is in / = 2, where the 
interaction is weakly repulsive. Thus, the effect of the 
imaginary part is considered to be small for the decay 
channels of Eqs. ©, ([3]), and Qj. This is indeed seen 
in the numerical estimation in Fig. [51 In any event, to 
construct a reliable model for the Ac — ?> ttttE decay, de- 
tailed analysis of the experimental decay process will be 
needed. 



Appendix B: RELATIVE SIGN OF THE WEAK 
INTERACTION 

Here we present the estimation of the relative sign of 
the weak vertices using the quark diagram method [s^ . 
The weak decay of the c quark in Ac occurs through W'^ 
exchange as 



In addition, it is necessary to create qq pair(s) from the 
vacuum to obtain the ttttE final state. 

Let us first consider the decay processes Ac — > 
7r"'"(7r~E"'") in detail. Four diagrams shown in Fig. [HI can 
contribute to this process. As shown in Appendix 1X1 the 
TT"*" has a large momentum. This is favored when the tt"*" 
is created by the weak decay of the c quark. In this view- 
point, the contribution from Fig.lHl^d) is considered to be 
small. In addition. Figs. [IJb) and M,c) can be obtained 
by Fierz rearrangement from Fig. ^a,) which introduces 
a suppression factor, as we will see below. 

The matrix element of the weak process [Fig.[9l[a)] can 
be written as 



M 



(a) 

7r+(7r-E+) 



^(7r+(7r-E+) \{uu)[ur^d][sT''c]\A, 



where F^ = 7/^(1^75), we insert the operator of the weak 
interaction, and the uu pair is created from the vacuum. 
The primary pion 7r+ has large momentum so it can be 
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(a) 




(b) 

d ■ 



d' 
d\ 

u 2° 




(d) 




FIG. 9. Quark diagrams for the Ac — > ■k'^{-k E""") decay. 



factorized as 



M 



(a) 



-(7r+ |ur^d|0)(7r"I]+ |(uu)[sr'"c]| Ac 
--tkf U{7r-J:+ \{uu)[sr^c]\A,), 



where k'^ is the momentum of tt^ and f-^ is the pion 
decay constant. The 7r+ is created by the axial vector 
component of the matrix element. Since the tt^ has a 
small momentum, we can apply the soft pion theorem to 
obtain 



M 



(a) 



^kf{j:+\[Q~,{uu)][sT'^c]\A, 



(Bl) 



Before addressing the other decay channels, we turn to 
Fig.EJb). This amplitude can be obtained by interchang- 
ing the u quarks in the amplitude [Fig.|9l^a)]. Noting that 
the relevant terms to create 7r+ are the pseudoscalar and 
axial vector combinations in the Fierz rearrangement Q 



d' 
d 



^ I n'^ (high momentum) 



FIG. 10. Dominant quark diagram for the Ac — >■ 7r^(7r^E ) 
decay. 



* Here we assume that uu creation takes place at the same point 
with the weak interaction vertex, which is not always the case 
in reality. However, the suppression factor by color and spinor 
indices as well as the sign of the amplitude are not modified by 
the nonlocality of the operator. 



FIG. 11. Dominant quark diagrams for the Ac — s- 7r"'"(7r''E'') 
decay. 



we obtain 

(tt+tt" |(?2iu)(u2r,jd/j) 
==i(7r+7r" ||i(M27i'75rf)(wir^7''75u) 
- i('"275d)(uir^75M)|, 

where the suppression factor 1/3 is responsible to obtain 
the color singlet pair. Thus, the first term is suppressed 
by the factor 1/12 in comparison with the diagram (a). 
The second term creates the tt"*" from the pseudoscalar 
operator. The matrix element can be related with that 
of the axial vector current as 



!7r+ |u7^75d|0) 



= -(7r|il75d|0). 
X 



>+ 1^75^1 0) 



In the case of the hyperon nonleptonic decay, the en- 
hancement factor X is large enough to cancel the sup- 
pression factor by the Fierz rearrangement [s^ . In the 
present case, however, kj^^ is the momentum of the pri- 
mary pion {\kf*\ 747 MeV and fcj^ - 760 MeV) so 
the enhancement factor is 



X 



mz 



2.6, 



and the pseudoscalar component is not much enhanced. 
Thus, we consider Fig. [Sl^a), in which 7r+ is created from 
the decay of the boson, provides the main contribu- 
tion to this weak process. 

Next, we consider the Ac 7r+(7r+S~) process. The 
diagram corresponding to Fig. is shown in Fig. [TU] 
where the tt"*" generated in the decay is regarded as 
the primary pion. In this case, the dd pair should be 
created from the vacuum instead of uu. Following the 
same step as before, we obtain 



M 



7r+{7r+i;-) 



;s]- |J75M[sr'^c]|Ac 



(B2) 



The dominant contributions to the Ac — > 7r+(7r S ) pro- 
cess are shown in Fig. [TT] We have two analogous dia- 
grams with uu creation [Fig. Illf a)] and with dd creation 
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TT (high momentum) 



FIG. 12. Quark diagram for Ac 7r"(7r+E°) decay. 



_ I 7r (high momentum) 

d 



FIG. 13. Quark diagram for Ac -> 7r"(7r"E+) decay. 
[Fig. [TlT b)]. From these diagrams, we obtain 




(B3) 



For the processes Ac — )■ 7r°(7r+E°) and Ac -> 
TT*^ (71*^2+), in which the 7r° is the primary pion, there 
is no corresponding diagram to Figs. [10] and [TTJ Thus, 
the leading contribution would be one of the Fierz rear- 
ranged diagrams. Here we estimate the relative sign by 
taking the graphs in Figs. [T2land[T3l which can be written 
as 

X.o(,+so) ^(7r"(^+S°)|(dd)[Mr^d][sr^c]|Ac), 
X,o(,oso) ^(7r"(^0l]+) |(Jd)[Mr^d][sr^c]| Ac). 

Utilizing the Fierz transformation and soft pion theorem, 
we obtain 

^2 

■A^irO(7r+EO) ' 



S+) 



12(to„ + ma) 

(s°|[sr^u- Jr^rf][sr^c]|Ac), (B4) 
;s+|i2r^d[sr^c]|Ac), (B5) 



12 (m„ + md) 



where we have neglected the relatively small axial vector 
components. 

In this way, we obtain the main components of the ma- 
trix elements of the weak decay processes. To determine 
the relative sign, we rotate the matrix elements in the 
isospin space. The phase convention for the isospin state 
is given by 



\7T+) 



\7TiI=l,l3-- 



(B6) 



as well as a negative sign for dj^u ~ 7r+. For instance, it 
follows from Eq. ()Bip that 



M 



(a) 

7r+(7r-S+) 



k; (s+|s75d[sr'^c]|Ac) 

= -kf{^-\dj,u[sT''c\\A,), 



where the arrow indicates the isospin rotation. Compar- 
ing this with Eq. (|B2p we find that the amplitude of the 
Ac 7r+(7r^E+) process has the same sign with that 
of the Ac — > 7r+(7r+I]~) process. In the same way, from 
Eq. ((B3| we obtain 



X.+ (.oso) + (E-ld^sT/isF^cllAc). 

In comparison with Eq. (|B2p , we find the relative sign of 
Ac ^ 7r+(7r-I]+) and Ac ^ 7rO(7rOl]+) is odd. Finally, 
from Eq. (|B4[) we obtain 



M 



7rO(7r+S") 



12(m„ -I- md) 



;s+|[i2r^d][sr^c]|Ac 



When we compare it with Eq. (|B5[) , we find the same sign 
of the processes Ac 7r°(7r+E°) and Ac tt°{tt°T.+). 
In this way, we determine the relative sign of the weak 
interaction vertices in Table |lTl 



[1] Y. Nogami, Phys. Lett. 7, 288 (1963). 
[2] Y. Akaishi and T. Yamazaki, Phys. Rev. C 65, 044005 
(2002). 

[3] N. V. Shevchenko, A. Gal, and J. Mares, Phys. Rev. Lett. 

98, 082301 (2007). 
[4] Y. Ikeda and T. Sato, Phys. Rev. C 76, 035203 (2007). 
[5] N. V. Shevchenko, A. Gal, J. Mares, and J. Revai, Phys. 

Rev. C 76, 044004 (2007). 
[6] T. Yamazaki and Y. Akaishi, Phys. Rev. C 76, 045201 

(2007). 



[7] A. Dote, T. Hyodo, and W. Weise, Nucl. Phys. A 804, 
197 (2008). 

[8] A. Dote, T. Hyodo, and W. Weise, Phys. Rev. C 79, 
014003 (2009). 

[9] Y. Ikeda and T. Sato, Phys. Rev. C 79, 035201 (2009). 
[10] Y. Ikeda, H. Kamano, and T. Sato, Prog. Theor. Phys. 

124, 533 (2010). 
[11] T. Hiraiwa et al. (J-PARC E15 Collaboration), Int. J. 

Mod. Phys., A 26, 561 (2011). 
[12] K. Suzuki et al. (FOPI Collaboration), Nucl. Phys. A 



13 



827, 312C (2009). 
[13] J. Zmeskal, M. Bazzi, M. Bragadireanu, P. Buhler, 

M. Cargnelli et al, Nucl. Phys. A 835, 410 (2010). 
[14] M. Bazzi et al. ,arXiv:1105. 30901 [nucl-ex]. 
[15] N. Kaiser, P. B. Siegel, and W. Weise, Nucl. Phys. A 

594, 325 (1995). 
[16] E. Oset and A. Ramos, Nucl. Phys. A 635, 99 (1998). 
[17] J. A. Oiler and U. G. Meissner, Phys. Lett. B 500, 263 

(2001). 

[18] M. F. M. Lutz and E. E. Kolomoitscv, Nucl. Phys. A 

700, 193 (2002). 
[19] T. Hyodo and D. Jido, Prog. Par t. Nucl. Phys. 

(2011), doi:10.1016/j.ppnp. 2011. 07.02, larXmT0lL4474l 

[nucl-th] . 

[20] D. Jido, J. A. Oiler, E. Oset, A. Ramos, and U. G. Meiss- 
ner, Nucl. Phys. A 725, 181 (2003). 
[21] T. Hyodo and W. Weise, Phys. Rev. C 77, 035204 (2008). 
[22] B. Borasoy, R. Nissler, and W. Weise, Phys. Rev. Lett. 

94, 213401 (2005). 

[23] B. Borasoy, R. Nissler, and W. Weise, Eur. Phys. J. A, 
25, 79 (2005). 

[24] J. A. Oiler, J. Prades, and M. Verbeni, Phys. Rev. Lett. 

95, 172502 (2005). 

[25] J. A. Oiler, Eur. Phys. J. A, 28, 63 (2006). 

[26] B. Borasoy, U. G. Meissner, and R. Nissler, Phys. Rev. 

C 74, 055201 (2006). 
[27] Y. Ikeda, T. Hyodo, D. Jido, H. Kamano, T. Sato, and 

K. Yazaki, Prog. Theor. Phys., 125, 1205 (2011). 
[28] J. M. Blatt and V. F. Weisskopf, Theoretical Nuclear 

Physics (John Wiley & Sons, New York, 1952) 
[29] G. Beer et al. (DEAR Collaboration), Phys. Rev. Lett. 

94, 212302 (2005). 
[30] B. Adeva et al. (DIRAC Collaboration), Phys. Lett. B, 

619, 50 (2005). 
[31] J. Gasser, V. Lyubovitskij, and A. Rusetsky, Phys. Rep. 

456, 167 (2008). 
[32] J. R. Batley et al. (NA48/2 Collaboration), Eur. Phys. 

J. C 54, 411 (2008). 



[33] N. Cabibbo, Phys. Rev. Lett. 93, 121801 (2004). 
[34] N. Cabibbo and G. Isidori, J. High Energy Phys., 03 
(2005) 021. 

[35] U. G. Meissner, G. MuUer and S. Steininger, Phys. Lett. 

B 406, 154 (1997); 407, 454 (1997). 
[36] J. R. Batley et al. (NA48/2 Collaboration), Phys. Lett. 

B 633, 173 (2006). 
[37] J. R. Batley et al. (NA48/2 Collaboration), Phys. Lett. 

B 686, 101 (2010). 
[38] S. Barlag et al. (ACCMOR Collaboration), Phys. Lett. 

B 283, 465 (1992). 
[39] Y. Kubota et al. (CLEO Collaboration), Phys. Rev. Lett. 

71, 3255 (1993). 
[40] P. Avery et al. (CLEO Collaboration), Phys. Lett. B 325, 

257 (1994). 

[41] P. L. Frabetti et al. (E687 Collaboration), Phys. Lett. B 

328, 193 (1994). 
[42] E. Vazquez- Jauregui et al. (SELEX Collaboration), Phys. 

Lett. B, 666, 299 (2008). 
[43] K. Nakamura (Particle Data Group Collaboration), J. 

Phys. G 37, 075021 (2010). 
[44] Y. R. Liu and S. L. Zhu, Phys. Rev. D 75, 034003 (2007). 
[45] M. Mai, P. C. Bruns, B. Kubis and Ulf-G. Meissner, Phys. 

Rev. D 80, 094006 (2009). 
[46] T. Hyodo, D. Jido, and A. Hosaka, Phys. Rev. Lett. 97 

192002 (2006). 

[47] T. Hyodo, D. Jido, and A. Hosaka, Phys. Rev. D 75, 
034002 (2007). 

[48] A. Torok, S. R. Beane, W. Detmold, T. C. Luu, 
K. Orginos, et al, Phys. Rev. D 81, 074506 (2010). 

[49] J. A. Oiler, E. Oset and J. R. Pelaez, Phys. Rev. D 59, 
074001 (1999). 

[50] M. Fukugita, Y. Kuramashi, M. Okawa, H. Mino, and 
A. Ukawa, Phys. Rev. D 52, 3003 (1995). 

[51] Y. Ikeda (HAL QCD Collaboration), Prog. Theor. Phys. 
Suppl. 186, 228 (2010). 

[52] L. B. Okun, Leptons and Quarks (North-Holland, Ams- 
terdam, 1982). 



